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TECHNICAL MEMORANDUM 


THEORETICAL REGIME DIAGRAMS FOR THERMALLY DRIVEN 
ROWS IN A BETA-PLANE CHANNEL IN THE 
PRESENCE OF VARIA8LE GRAVITY 

1. INTRODUCTION 


It is generally agreed that the wave-like disturbances seen in 
rotating cylindrical geometry laboratory flow experiments arc the result 
of baroclinic instability of axially symmetric flow. This interpretation has 
been given a good theoretical basis by the study of Burcilon [1]. in which 
the Eady model of baroclinic instability was used to obtain stability criteria 
for axially symmetric flows. These criteria are in reasonable agreement 
with those observed in the laboratory for rotating annulus flows. 

in the Earth's atmosphere, baroclinic instability is also an important 
process for maintaining departures from axially symmetric flow. It has 
long been believed that better simulation of atmospheric flow patterns, or 
at least a better understanding of how and when baroclinic instability 
operates on the atmosphere, could be achieved if laboratory rotating fluid 
experiments could be done in spherical geometry. Such experiments have 
not been realizable because the dielectric body force for simulated radial 
gravity cannot be made large enough to dominate the effect of ambient 
terrestrial gravity in the laboratory. The low gravity environment aboard 
orbiting laboratories such as Spacelab, to be operational in the ^arly 
1980's, affords an opportunity for such an experiment. 

In going from cylindrical geometry to spherical geometry in a rotating 
fluid experiment, one important new feature is the latitudinal variation of 
the local vertical component of rotation. As is well known, the effect of 
this on the dynamics of low -frequency geophysical motions can be taken 
into account by ; ; -plane geometry. As one of the first steps in developing 
a model for use in the design of a Spacelab experiment, Geislcr and Fowlis 
{?] extended the work of Barcilo.i [1] to a t^-plane channel. The princi- 
pal result of their study was to document the changes in the shape and 
location of the baroclinically unstable region of parameter space brought 
about by the latitudinal dependence of the vertical component of rotation. 

One consequence of using a dielectric body force to simulate gravity 
is that the force field law is one of inverse fifth power [3]. This must 
be taken into account in mathematical models of the proposed experiment 
and, moreover, is potentially troublesome because it does not simulate the 
inverse square of terrestrial gravity. 



This report describes the extension of the barocUnic instability 
model of Geisler and FowUs 121 to include an inverse fifth power law of 
gravity. The study shows that there is little difference between the sta- 
bility information obtained from the two models, provided results are 
plotted using a vertical average of gravity in the plotting parameters. 
This result supports conclusions obtained earlier by Giere and FowUs (41. 


II. FORMULATION 


BarocUnic instability in the presence of constant gravity was treated 
in Reference 2. We obtained growth rates and eigenfunctions for unstable 
modes in both the Charney and Eady models of barocUnic instability with 
and without Ekman damping at the boundaries. The normal modes were 
assumed to have the functional form 


i|<x,y,z,t) = 4>(z) sin 



exp lik(x-ct)] 


( 1 ) 


where ^ is a stream function, h is the width of the channel, and n is a 
positive integer. The equation solved there was 


o/V d\j 6 _L 

a y + |(TFT> 


subject to the boundary conditions 


I ilf n v d dU , N 2 

ik (U " c) S ‘ 1 W 

\ I J 0 


where the plus sign appUes at the upper boundary and the minus sign 
applies at the lower boundary. 

In the preceding equations, U(z) is the basic state flow whose 
stabiUty is being examined. As in Reference 2, we assume that the part 
of basic state temperature Geld T(y,z) associated with U(z) decreases 
Unearly with y and the basic state temperature field <T(z)> associated 
with the static stability of the fluid increases Unearly with z. These 
parameters enter the problem through the thermal wind equation 



2 



( 4 ) 



find through the definition of N 2 


N * * «« (tF) • < 5) 

Here T is the zonally averaged temperature, <T> is the area averaged 
temperature, a is the coefficient of thermal expansion of the fluid, and g 
is gravity. In this report we assign to g the variation 


g = 


*o 

(1 + z/a) p 


( 6 ) 


where p is an integer, a is the inner radius of the laboratory device, and 
g 0 is the value of g at z = 0. Integration of (4) upward from z = 0 

(where we trite U = 0) gives the basic state flow 


U(2 > - - -r- If 

o 


1 


(p - D 


1 - (1 + z/a) 


-p+1 


(7) 


If p = 1, integration of (4) gives the logarithmic flow 


U(z) = 


a *o a 3T r n 

W L (1 + z/a) ] 


(8) 


The stability of the flow given by equation (8) was examined by Giere 
and Fowlis [4]. In the present report we examine the stability of the 

flow given by equation (7) with p = 5, that is, an inverse fifth power 

gravity . 

In the case of U(z) more general than the linear variation with z- 

used in Geisler and Fowlis [2], the parameter fl in equation (2) should 

be replaced by 


3 



g d/ f o 2 du\ 

8 ' *7 ’ 


( 9 ) 


However, as can be seen from equations (4) and (5), the factor g cancels 
out and the correction to 6 then vanishes when 9T/3y and d<T>/dz are 
constant, as is the case here. 

We introduce nondimensional quantities denoted by a prime as 
follows : 


x/L 

k' = kL 



y/L 

U' = U/AU 


(10) 

z/d 

C* = C/AU 

a 



Here L is an arbitrary horizontal length scale taken to be 0.707a in 
Reference 2, and d is the depth of the fluid. The quantity AU is U(d), 
that is, the difference in the basic state flow between z = 0 and z s d. 
Equations (2) and (3) then become 



The parameters in these equations are a static stability S, a 6-parameter 
B , and Ekman number E and a thermal Rossby number R . They are 
defined as ° 


S 


N 2 d 2 



(13) 


E 




(14) 
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In the sequel, we refer to models with B * 0 in equation (12) as 
inviscid models. We refer to equations (11) and (12) as the Charney 
model of baroclinic instability. We refer to the model obtained by setting 
B = 0 in equation (11) as the Bady model. 

In the inviscid Bady model it is customary to display the behavior 
of Cj (the imaginary part of the mode phase speed) in a diagram of 

KS^ 2 Cj versus KS 1 ^ 2 . Here S is given by equation (13) and 



Giere and Fowlis [4] studied the unstable modes in the inviscid Eady model 
using the log profile basic state flow given by equation (8). For the value 
of AU in equation (10) that nondimensionalizes U and c, they used the 

— T73 2 

quantity obtained by replacing g with (g in equation (4) and then 
integrating this equation from 0 to d. They found that the curve changed 

very little over a wide range of the parameter d/a if is used in the 

- 1/2 T77 

plot instead of S , where S is the value of S obtained from using 
in the definition of N [see equations (6) and (13)] the vertical averaged 
1/2 

value of g ' defined as 




(16) 


Unless otherwise stated, in the present study we have used 

AU = U(d) to nondimensionalize U and c, and we have used instead of 
777 _i/o 

S the quantity S , where S is the value of S defined using the verti- 

l /2 

cal average of g rather than g . 


III. INVISCID MODEL RESULTS 


All results presented here imd in Section IV wore obtained by solv- 
ing equation (11) subject to boundary conditions in equation (12) using 
the numerical technique briefly described by (leislcr and Fowlis [ 2) . As 
noted in Section II , the difference between that and the present study is 
that we here use gravity which varies with radial distance according to 
the function 


5 



where a is the inner radius and d the depth of the fluid. 

We consider first the case Of p ■ 1 in the in viscid Bady model, the 
case for which analytic solutions were obtained by Giere and Fowtts [4]. 

We show in Figure i a plot of KS 1 ^ 2 c'j for the runs with d/a = 0 (that 

is, oonstant g), d/a = 1, and d/a = 10. The figure shows that inverse 
first power gravity has little effect when the vertical average gravity is 

1/2 

used for the plot. When the vertical average of g is used, as by 
Giere and Fowlis [4], the three curves for d/a * 0, 1 and 10 cannot be 
distinguished from one another. This agreement with the analytic results 
of Giere and Fowlis [4] indicates that the numerical routine is functioning 

I/O 

correctly and also confirms that the vertical average of g rather than 
g is optimum for that problem. 

Figure 2 shows results for the inviscid Eady model when p = 5 
(inverse fifth power gravity). Runs shown are those for d/a = 0, 
d/a = 0.2 and d/a = 1. We have not gone beyond d/a = 1 because we 
anticipate that d/a will have to be <1 in any reasonable geophysical experi- 
ment. Figure 3 shows results when we go to the in viscid Charney model 
(adopting the value B = 2.35 used as a standard case by Geisler and 
Fowlis [2]). Figures 2 and 3 support the conclusion that, for the range 
of d/a of geophysical interest at least, such a plot is little affected by 
inverse fifth power gravity if the vertical average gravity is used. If 

1/2 

the vertical average of g is used (together with the corresponding 
change of definition of AU), the position of the curves relative to the 
d/a = 0 case changes slightly. The results for the inviscid Eady model 
and the inviscid Charney model are shown, respectively, in Figures 4 

i/4 —1/2 

and 5. These show that the use of S rather than S is still pref- 
erable in the Eady model, but not in the Charney model. 


IV. REGIME DIAGRAMS 


Superimposed contour diagrams of c. for several values of sonal 

* 2 2 

wave number k in the parameter space of S versus R Q / S E constitute a 

theoretical regime diagram . The construction of these diagrams for the 
Eady model and the Cnarney model was the subject of the paper by Geisler 
and Fowlis [2]. Figure 6 is taken from that paper and shows a theoret- 
ical regime diagram for the Eady model when d = h * L * 0.707a. 


(The parameter e 2 / 5 = R q 2 /S 2 E. ) The curve labelled (1) is the stability 

boundary for zonal wave number (1) disturbances, outside (to the left) 
of which Cj < 0. The envelope of the curves shown in this figure 

obviously separates the region of parameter space where the flow is 
unstable from that where the flow is stable and, hence, axially symmetric. 
Geophysical experiments must be such that the unstable regime can exist 
in the apparatus, hence the utility of theoretical regime diagrams in 
experimental design studies. 

To illustrate the results of our study of the effect of inverse fifth 
power gravity on theoretical regime diagrams, we have selected zonal 

- 1/2 

wave number 3. Figure 7 shows contours of KS e'j for this wave num- 
ber in the Eady model when g is constant. The second case, where 
d/a = 0.2, is so close to this that it does not bear showing. Figure 8 
shows the extreme (for a geophysical experiment) case for d/a = 1. 
Comparison with Figure 7 shows that the difference is rather small. We 
have run cases for other zonal wave numbers with the same result. The 
conclusion is that for d/a < 1, the inverse fifth power gravity has only 
a small effect on the shape and location of the unstable regime provided 
vertically averaged g is used in drawing the diagram. 

- 1/2 

Figure 9 shows contours of KS ' c'j for zonal wave number 3 in the 

Charney mode 1 when g is constant. Here we adopt the value B = 2.35 
used as a standard case by Geisler and Fowlis [2]. The case for 
d/a = 0.2 is again not much different and is not shown here. Figure 10 
shows the case for d/a = 1. The change is somewhat greater than in the 
corresponding Eady model runs, the most notable change occurring in the 
region S < 0.1. However, there is very little change in the leftward 
penetration of the nose located at about S = 0.2. We have obtained similar 
results for other* zonal wave numbers. We conclude that for d/a < 1 the 
inverse fifth power gravity does not have significant effect on regime 
diagrams in either the Charney model or the Eady model provided that 
parameters based on a vertically averaged gravity are used. 
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Figure 1. Imaginary part of phase speed in the in viscid Eady model 
for the case of inverse first power gravity. Curves are labelled 
by the value of d/a. The dashed curve is the case of 
constant gravity. 



Figure 2. Imaginary part of phase speed in the in viscid Eady model 
for the case of inverse fifth power gravity. Curves are labelled by 
the value of d/a. The dashed curve is the case of 
constant gravity. 
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Figure 3. Imaginary part of phase speed in the inviscid Charney 
model for the case of inverse fifth power gravity. Curves 
are labelled by the value of d/a. The dashed curve is 
the case of constant gravity. 
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Figure 4. Same as Figure 2, but with S ' as parameter. 
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Figure 5. Same as Figure 3, but with S as parameter. 
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Figure 6. Theoretical regime diagram for the Eady model with 
d = h = accs0 o (from Geisler and Fowlis [2]). The curves 

are stability boundaries labelled by zonal wave number. 
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Figure 7. Contours of imaginary part of phase speed KS c^' 
in the Eady model with gravity constant. 



- 1/2 

Figure 8. Contours of imaginary part of phase speed KS Cj* in 

the Eady model with inverse fifth power gravity and 

d/a = 1. 
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